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' In this paper we investigate the influence of an external magnetic field on a flavoured 

holographic gauge theory dual to the D3/D5 intersection at finite temperature. Our 
study shows that the external magnetic field has a freezing effect on the confinement/ 
deconfinement phase transition. We construct the corresponding phase diagram. We 
investigate some thermodynamic quantities of the theory. A study of the entropy reveals 

enhanced relative jump of the entropy at the "chiral" phase transition. A study of the 



> 

l/^ magnetization shows that both the confined and deconfined phases exhibit diamagnetic 

jn response. The diamagnetic response in the deconfined phase has a stronger temperature 

dependence reflecting the temperature dependence of the conductivity. We study the 
meson spectrum of the theory and analyze the stability of the different phases looking 
at both normal and quasi-normal semi-classical excitations. For the symmetry breaking 
O phase we analyze the corresponding pseudo-Goldstone modes and prove that they satisfy 

J> non-relativistic dispersion relation. 

X 
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1 Introduction 



Recently a lot of attention has been focused on the properties of flavoured holographic gauge 
theories and their potential application to realistic phenomenological models. One of the main 
motivations for such studies is the simple geometric description that holographic models provide 
of non-perturbative phenomena difficult to study via conventional field theoretical methods. 

Despite the growing spectrum of applications of the AdS/CFT correspondence [T] a major 
limitation is that realistic field theories do not seem to have simple holographic backgrounds. 
This is why it is of particular importance to investigate properties of non-abelian gauge theo- 
ries exhibiting universal behaviour. Particularly interesting is to analyze the phase structure 
of strongly coupled Yang-Mills theories. An example of such application of the holographic 
approach is the study of properties of strongly coupled quark-gluon plasmas (see refs. O ]3] for 
a recent review). 

The rich and complicated dynamics of strongly coupled non-abelian gauge theories suggests 
investigation of their response to strong controlling parameters such as temperature, external 
electromagnetic fields and various chemical potentials. In this regime one may expect that 
otherwise different theories would exhibit similar behaviour and one can extract valuable results 
by analyzing holographic gauge theories. This is the context of our present studies. 

Our holographic set up employs holographic gauge theories dual to the Dp/Dq brane intersec- 
tion. This class of gauge theories has been extensively studied in the literature. The first and 
most studied example is the D3/D7 intersection [1]. Early studies of the meson spectrum have 
been performed in ref. 0. The phase structure of the theory at finite temperature has been 
studied in [H] and a confinement / deconfinement phase transition corresponding to topology 
changing transition of the probe D7-brane has been revealed. More detailed studies of this 
phase transition have been performed in refs. [7]-[T3j. The phase transition has been classified 
as a first order one. Universal behaviour of the general Dp/Dq system has been uncovered. 
Various properties of the general Dp/Dq set up in constant external electric and magnetic fields 
have been analyzed in refs. [20]-|34]. 

One major aspect of the analysis presented in this paper is the phenomena of mass generation 
in an external magnetic field. This phenomena has been extensively studied in the conventional 
field theory literature [Il]-[l9]. The effect was shown to be model independent and therefore 
insensitive to the microscopic physics underlying the low energy effective theory. The essence 
of this effect is the dimensional reduction D D — 2 (3 + 1 — > 1 + 1) in the dynamics of fermion 
pairing in a magnetic field. 
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Magnetic catalysis of mass generation has been demonstrated in various 1+2 and 1+3 dimen- 
sional field theories. Given the universal nature of this effect it is natural to explore this 
phenomena in the context of holographic gauge theories. Such studies for the Dp/Dq inter- 
sections have been extensively performed for both 1+3 and 1+2 dimensional theories. In ref. 
j31j an universal description of the general Dp/Dq system has been attempted. It has been 
demonstrated that the mechanism of spontaneous symmetry breaking in external magnetic field 
for the Dp/Dq system is a universal feature of this class of theories. 

Particularly interesting are the findings for the D3/D5 set up. Recently this theory received a 
great deal of attention and emphasis has been made of the potential application of this brane 
configuration in describing qualitative properties of 1 + 2 dimensional condensed matter systems 
(see for example refs. [29], [30l |35]). In ref. [31] it has been shown that the pseudo goldstone 
modes of this theory satisfy non-relativistic dispersion relation. The low energy effective action 
has been obtained. To the computed order, it agrees to the effective action of spin waves in a 
ferromagnetic. The existence of a single time derivative term in the effective action confirms 
the potential phenomenological applications of this model. 

In this paper we extend our previous investigation of the D3/D5 set up in an external magnetic 
field to the case of finite temperature. However, we are not only interested in expanding the 
description of the pseudo-goldstone modes to the finite temperature case but also to study the 
effect that the external magnetic field has on the confinement /deconfinement phase transition. 
This is why we are interested in both the small and large bare mass sectors of the theory. 
Similar studies of the D3/D7 set up have been performed in [Ml [25] . It has been shown 
that the magnetic field has a freezing effect on the phase transition and for sufficiently strong 
magnetic field the deconfined phase cease existence (for any bare mass). Our studies reveal the 
same qualitative behaviour in the case of the D3/D5 set up. However there are some differences 
that are of potential interest. 

Let us summarize the content of the paper: 

Section 2 of this work is separated into two subsections. In the first subsection we provide a brief 
description of the holographic set up. The parameterization of the background geometry that we 
use, the basic extracts from the AdS / CFT dictionary needed for the holographic description as 
well as the method of introducing external magnetic field, are presented. The second subsection 
studies the equation of state of the system in the condensate versus bare mass plane. The study 
involves both analytic and numerical techniques. The analytic results are for large bare masses. 
An asymptotic expression for the condensate as a function of the magnetic field, the temperature 
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and the bare mass is obtained. It is demonstrated that unhke the 1+3 dimensional case (the 
D3/D7 set up) in the 1+2 case the competing effect of the magnetic field and the temperature 
is present even at large bare masses. This gives the possibility to tune the condensate of 
the theory to zero for arbitrarily large masses. The reason for that could be that unlike in the 
D3/D7 set up the external magnetic field does not require new counter terms in the holographic 
regularization of the D-brane action. For small bare masses we perform a numerical analysis. 
The freezing effect of the magnetic field is demonstrated. The existence of critical ratio of 
the temperature and the square root of the magnetic field beyond which the phase transition 
disappears is demonstrated. 

Section 3 of this parer is dedicated to the thermodynamic analysis of the theory. A holographic 
renormalization of the D5-brane action in the spirit of p7] is performed and the free energy 
of the theory is obtained. The proper thermodynamic ensemble is identified and the phase 
diagram of the theory in various coordinates is constructed. The temperature dependence of 
the condensate at zero bare mass is explored as an order parameter for the "chiral" phase 
transition. A study of the entropy is performed. It is demonstrated that at zero bare mass the 
jump of the entropy is enhanced relative to the jump of the entropy at large bare masses. A 
possible explanation is that at zero bare mass the confinement /deconfinement phase transition 
is also a chiral phase transition. Finally, the magnetization of the theory is computed. It is 
demonstrated that both the confined and deconfined phases are diamagnetic. In the confined 
phase the diamagnetic response of the theory is almost independent on the temperature. At the 
phase transition the magnetization has a large negative jump and a very strong diamagnetic 
response possibly because the theory is in a conductive phase. Interestingly the diamagnetic 
response of the deconfined phase depends strongly on the temperature and vanishes inversely 
with the temperature. The strong temperature dependence reflects the temperature dependence 
of the conductivity. At large temperature the susceptibility is deduced by dimensional analysis 
reflecting the return to conformality. 

Section 4 of this work studies the meson spectrum of the theory. First the spectrum of fluc- 
tuations corresponding to the radial coordinate in the transverse are considered and the 
stability of the theory is analyzed. Both normal and quasi-normal excitations are studied. It 
is demonstrated that the thermodynamically stable phases are tachyon free. It is also verified 
that phases with negative heat capacity have tachyonic spectrum and are thus unstable. The 
existence of metastable phases accessible by supercooling is inferred. At zero bare mass the 
case of zero magnetic field and finite temperature is analyzed along the lines of ref. [12]. The 
corresponding Heun equation is solved employing the method of continued fractions. Next, the 
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spectrum corresponding to fluctuations along the angular coordinates of the transverse is 
analyzed. Only the confined phase is considered. The qualitative features of the spectrum are 
the same as for the zero temperature case studied in pIT]. At large bare mass a Zeeman splitting 
of the energy levels is observed. For small bare mass the pseduo-Golstone mode is analyzed. A 
non-relativistic dispersion relation is observed. The analysis of the pseudo-Goldstine modes is 
performed both numerically and analytically. 

We end the paper with a brief conclusion in Section 5. 



2 Preliminaries 
2.1 General set up 

Let us consider the supergravity background corresponding to the near horizon limit of Nc 
coincident D3-branes at finite temperature: 

_ 1 ^(4) _ 1 (4r^ + 

where (IVL2 = da^ + cosa^d/^^; dVL\ = dip^ + cos^ = + f'; . 

The metric in equation ([T]) corresponds to the AdSs x black hole in a Poincare patch. Its 
extremal limit is holographically dual to an A/" = 4 SU{Nc) Supersymetric Yang-Mills theory 
in 1 + 3 dimensions [1]. The radius of the asymptotically AdSs space is related to the t'Hooft 
coupling of the gauge theory via = 2Xa'. The Hawking temperature of the AdSs black hole 
is given by T = h/nB^ and is interpreted as the temperature of the holographically dual gauge 
theory. For the choice of radial coordinate considered in equation ([T| the subspace transverse 
to the volume of the background colour branes is conformally M^. We have further split the 
transverse subspace to the product M? x and introduced spherical coordinates p, VL2 and 
1,0.2 in the corresponding subspaces. 

Next we introduce Nf probe D5-branes extended along the t, xi, 0:2, p, directions of the 
geometry. In the extremal limit such embeddings preserve half of the original supersymmetries 
of the background. The lowest energy sector of strings stretched between the D3 and D5 
branes gives rise to Nj fundamental M = 2 hypermultiplets confined on a 1 + 2 dimensional 
defect [HI El Hn]- The asymptotic separation of the D3 and D5 -branes in the transverse M? 
subspace is parameterized by / and is proptional to the bare mass of the fundamental Af = 2 
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hypermultiplets. If we consider the following anzatz for the D5-brane emebedings: 

/ = l{p) ; V = ; = , (2) 

the exact relation is given by: 

m = l{oo) = {2Tia)mq , (3) 

where is the bare mass of the fundamental field in the corresponding gauge theory. If the 
D3 and D5 branes overlap, the corresponding fundamental fields are massless and the theory 
has a global SO{?))r x 50(3) symmetry corresponding to independent rotations in the two 
subspaces the transverse space. A non-trivial profile of the D5-brane (Z(r) ^ 0) in 
the transverse M? subspace would break the global symmetry down to S'0(3)/j x f/(l), where 
U{1) is the little group in the transverse M? subspace. Furthermore, if the D5-brane wraps 
a shrinking S"^ cycle and closes at some radial (coordinate) distance above the horizon this 
distance is interpreted as the dynamically generated mass of the theory [13]. 

However the gauge theory dual to the D3/D5 system is conformal at zero bare mass and 
does not exhibit spontaneous symmetry breaking. One way to catalyze spontaneous symmetry 
breaking is to introduce a constant magnetic field to the theory. The effect of mass generation 
in external magnetic field has been widely studied on field theory side refs. [ll]-[19]. In the 
context of AdS/CFT the effect of magnetic catalysis in holographic gauge theories has been 
studied in [2n]-[S2]- In order to introduce magnetic field to the holographic set up we turn on 
a constant 5- field along the Xi,X2 directions of the geometry: 

B{2) = Hdxi A dx2 . (4) 

On field theory side this corresponds to introducing constant magnetic field H/ {27ra') per- 
pendicular to the plane of the defect. The D5-brane embedding is determined by the DBI 
action: 

Sbbi = -Nffi5 j rf'ee'*[-det(G,6 + B,, + 27ra'F,b)]V2 . (5) 

Ma 

Where Gab and Bab are the pull-back of the metric and the i?-field respectively and Fab is the 
gauge field on the D5-brane. 

With the anzatz ([2]) the Lagrangian is given by: 
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For large p':^ h the Lagrangian ^ asymptotes to: 



cos a , 



(7) 



which has a general solution of the form: 



l(p) = m H h . . . 

P 



(8) 



as we mentioned above the parameter m is the asymptotic separation of the D3 and D5 branes 
and corresponds to the bare mass of the dual gauge theory. The AdS / CFT dictionary suggests 
[1] that the parameter c is proportional to the condensate of the fundamental fields (ipip)- The 
exact relation is given by [7]: 



By solving the equations of motion of the probe D5-brane and extracting the asymptotic be- 
haviour at infinity one can generate the equation of state of theory in the condensate versus 
bare mass plane. This is one of the main tools that the holographic set up is providing and will 
be one of the main technique that we use to construct the phase diagram of the theory. 

2.2 The condensate of the theory 

In this section we study the fundamental condensate of the theory. We focus on the dependence 
of the condensate as a function of the bare mass for fixed temperature and magnetic field. We 
explore the effect of the magnetic field on the confinement /deconfinement phase transition of the 
fundamental matter. We show that the magnetic field has a freezing effect on the transition. We 
also show that for sufficiently strong magnetic field the theory develops a negative condensate 
at zero bare mass and hence the global SO (3) symmetry is spontaneously broken. We will first 
focus on the regime of large bare masses m ^ b, R\/H, which can be treated analytically. 

2.2.1 Exact results for large bare masses. 

As we discussed above large bare mass correspond to large separation between the color and 
flavour branes. In this regime the probes are nearly BPS objects, their embeddings are very 
close to the trivial / = one and we can expand: 




(9) 



(10) 



where ^(p) ^ m and a 



l{0) so that we have: 



m = e'(o) = . 



(11) 
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On the other side the radial coordinate r ^ b remains large along the D5-brane embedding 
hence we can expand the lagrangian ^ in powers of b/r. To leading order we obtain the 
following equation of motion: 

O.i^p)) - ap^ [^l^^ + ^l^T^. J + • • • (12) 



After solving for ^{p), imposing the boundary conditions (11) and substituting in (12) we can 
extract the following expressions for the bare mass m and fundamental condensate c: 

" = "+^^+0^,); (13) 



7i{b^ - AH^R^) 5n{5b^ + m^H^R"^ + IQH^R^) / 1 
32^^2 4096a6 ^ 1^ 



The first equation in (13) can be easily inverted to obtain our final expression for the funda- 
mental condensate as a function of the bare mass valid for large m: 

Tiib^-AH^R^) 7r(57b^ -56b^H^R^ + 592H^R^) I \ 

^ = ,2m^ + ^ 40%;;^^ ^""[m^)-^ ^''^ 

A few comments are in order: 

It is interesting that even for large bare masses the competition between the effect of magnetic 
field and finite temperature is still present and the condensate can be either negative or positive 
depending on the ratio b/R^/H. This is to be compared to the 1 + 3 dimensional case, studied in 
PHEH] (using the D3/D7 set up), where the leading contribution to the condensate is negative 
and proportional to R^H'^/m. As we are going to see later this seems to be related to the fact 
that in the 1 + 2 dimensional case the external magnetic field does not lead to new divergences 
in the on-shell action, while in the 1 + 3 dimensional case there is a logarithmic divergence due 
to the magnetic field |25j . 

To extend our analysis to the regime of large temperature and magnetic field we need to employ 
numerical techniques. 



2.3 Numerical results 

In this section we solve numerically the equations of motion for the D5-brane embedding derived 
from ^ and explore the dependence of the condensate on the bare mass for fixed ratio of the 
external magnetic field and temperature of the system. We show that qualitatively the behavior 
of the system is the same as in 1 + 3 dimensions for the D3/D7 setup [211 125] . 
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Note that the background geometry ([T]) has a horizon at Tq = 6^/2 and the probe D5-brane 
embeddings spht into two classes: The first class are Minkowski embeddings which wrap a 
shrinking 5*^ in the part of the geometry and close at some finite distance above the balck 
hole. The second class of embeddings are Black hole type of embeddings that fall into the black 
hole and have an induced horizon on their worldvolume. The transition from one class of em- 
beddings to the other class of embeddings is a topology change transition that corresponds to 
a first order confinement /deconfinement phase transition for the fundamental matter j6l [Tf [TT]. 
There is also a critical embedding separating the two classes of embeddings that has a conical 
singularity at the horizon of the geometry. A more detailed analysis of the dual gauge theory 
near the state corresponding to the critical embedding reveals a discrete self-similar behaviour 
of the theory. In this regime the theory is characterized by the existence of multiple (thermo- 
dynamically unstable) phases manifested by a double logarithmic structure of the equation of 
state curve in the condensate versus bare mass plane. We shall not analyze this unstable 
regime of the theory any further. Instead we shall focus on the influence of the external mag- 
netic field on the confinement /deconfinement phase transition. As we are going to show the 
effect of the magnetic field is to decrease the critical temperature at which this transition takes 
place. Furthermore we shall demonstrate that for sufficiently strong magnetic field the phase 
transition cease existence and Minkowski embeddings are the only thermodynamically stable 
ones. 

For our numerical analysis it is convenient to define the following dimensionless variables: 

p = p/m; ri = HR^/b^] m = m/b; c = c/b'^; l{p) = -l{p/b) = m + c/ p + . . . (15) 
Upon substitution into the Lagrangian ^ we obtain: 

1/2 / 1 a^2~4 \ 1/2 



1 \ / 1 \ / 1 IQf] r 



To solve numerically the equation of motion obtained from (|16j) we impose the following bound- 
ary conditions: 

1(0) = to; I'iO) = 0; for Minkowski embeddings, (17) 



Kp)\c.h. = io; Hp)\c.h. 

p 



; for Black Hole embeddings. 

c.h. 

Using equation (|8j) we can extract the condensate and bare mass parameters m, c. The resulting 
plots for 7] = 0,0.75,2,4.25 are presented in figure [T} One can clearly see a first order phase 
transition pattern. It is also evident that the critical bare mass m^r decreases as the magnetic 
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0.35 




Figure 1: As one can see the critical bare mass decreases as the magnetic field increases and for 
T]cr ~ 4.25 it vanishes. Beyond this point the lowest positive branch is the stable one and at ttt, = 
there is a negative condensate breaking the global SO{3) symmetry. 

field increases and for rjcr ~ 4.25 it vanishes. Beyond this point (for rj > rjcr) the lowest positive 
branch of the condensate vs. bare mass plot is the stable one. Furthermore at m = the theory 
develops a non-negative condensate and hence the global S0{3) symmetry is spontaneously 
broken. It turns out that this symmetry breaking phase exists even for smaller values of the 
magnetic field (3 < ?7 < 4.25). The analysis of the meson spectrum in Section 4 will show that 
this is in fact a metastable phase of the theory that becomes stable ai rj = rjcr- Note also that 
there are more than one phases with broken global symmetry and zero bare mass. However 
analysis of the free energy shows (see Section 3) that the phase corresponding to the lowest 
negative value of the fundamental condensate is the stable phase. 

For large ratios of the external magnetic field and temperature of the system [i] > rjcr) the 
deconfined phase is parametrically suppressed and the qualitative behaviour of the theory 
approaches the one studied in ref. [3l]. In this limit the condensate curve has a spiral structure 
and the stable phase of the theory corresponds to the lowest positive branch of the spiral. 

The critical mass fhcr in figure [T] was determined by evaluating the free energy of the system 
(alternatively we could use the equal areas law [25]), which suggests renormalizing the euclidean 
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on-shell action of the probe brane and analyzing the thermodynamic properties of the theory. 
Analyzing these properties and constructing the phase diagram of the theory will be the subject 
of the next section. 



3 Thermodynamic analysis 

Let us begin by calculating the free energy of the theory. 
3.1 Free energy 

Since our theory is at fixed temperature T and magnetic field B the density of of the thermo- 
dynamic potential describing the ensemble satisfies |25] : 

dF = -SdT - ndB . (18) 

Here 5* is the entropy density of the system, B = H/ (27ca') is the magnetic field and /x is the 
magnetization. 

Following ref. [13] we will identify the regularized wick rotated on-shell lagrangian of the D5- 
brane with the free energy of the theory. Let us introduce a cut-off at infinity, pmax, The wick 
rotated on-shell lagrangian is given by: 

^T" / 1 \ / 1 \ 1/2 / IP 2-4 \ 1/2 



/..^iV.^Ws^^y .py^l-^J^l+^J [l^J^.] V1 + /-, (19) 



mm 



where V3 = J d^x and l{p) is the solution of the equation of motion derived from (16). The 
parameter pmin = for Minkowski embeddings and pram = p|e.h. for Black hole embeddings. It 
is easy to check that for large pmax we have: 



^3 



(20) 

Note that the the constant S-field does not introduce new divergencies to the on-shell action. 
This is different from the 1 + 3 dimensional case studied via the D3/D7 set up, where a 
logarithmic divergency is introduced proportional to R^H^. The pj^ax divergency is present 
even at zero temperature and magnetic field and can be removed by the addition of appropriate 
boundary terms. These terms have been computed in [37], where an elegant renormalization 
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procedure for general Dp/Dq intersections was developed. One can easily check that in our 
coordinates the counter terms action boils down to — |Pmax- is then convenient to define the 
following finite dimensionless quantity: 



'D5 




1.3 



(21) 



Note that Jd5 implicitly depends on m trough the on-shell form of /(p). Our final expression 
for the density of the free energy is thenQ 

F = Nf^4nb'ij,,{m,r]') . (22) 
9s 



Using equation (22) we can calculate the free energy density of the dual gauge theory in the 
functional form F(T, _B,mg) = T^f{^,i^), therefore the quantity Ids oc F/T^ can be used 
to determine the preferred phase at a given temperature. It is straightforward to calculate 
the critical bare mass parameter fhcr at which the confinement /deconfinement phase transition 
takes place. Plots of /d5 vs. m for 77 = 0, 0.75, 3, 4.25 (the values from figure [l]) are presented in 
figure |2} The critical values rhcr are determined by the kinks of the free energy. If we refine our 
numerical analysis and scan through all possible values of 77 for witch there is a phase transition, 
namely < 77 < 4.25, we can generate the phase diagram of the theory. 



3.2 Phase diagram 

In order to construct the phase diagram of the theory we generate numerically the critical 
curve fhcr{f]) in the f] vs. m plane. The resulting phase diagram is presented in figure [s] The 
finite enclosed area corresponds to Black hole embeddings that have a qusi-normal mode ex- 
citations and hence correspond to the deconfined phase of the theory. The rest of the phase 
space correspond to Minkowski embeddings which are characterized by discrete normal exci- 
tations that are interpreted as meson-like bound states of the dual gauge theory. This is the 
confined phase of the theory. Across the critical curve the theory undergoes a first order con- 
finement /deconfinement phase transition. Note that for > ?7cr ~ 4.25 the phase transition 
disappears and the theory is in the confined phase. If we consider the horizontal axes of the 
diagram (m = 0), the confinement /deconfinement phase transition is a "chiral" phase transi- 
tion. For rj < rjcr the theory has an SO (3) global symmetry and no mass-gap, while for rj > 

^Note that to avoid repetition we will sometimes refer to the densities of the thermodynamic quantities of 
interest by omitting the term density 
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71=0.75 




Figure 2: Plots of Ids vs. m for 77 = 0, .75, 3, 4.25. The critical value fhcr is determined by the position 
of the kink of Ids • 



the theory is in a spontaneously broken phase with finite mass-gap. There are also Goldstone 
modes associated to the spontaneous symmetry breaking which we shall study in details in 
Section 4 when we focus on the meson spectrum of the theory. Qualitatively the phase diagram 
of the theory is the same as the one for the 1 + 3 dimensional case studied in [211 12^ . 

To facilitate comparison to phase diagrams obtained via non-holographic techniques it is useful 
to represent our phase diagram in field theory units. To this end note that the definitions of 
our dimensionless quantities m, 1] imply the following relations: 



1 

771 



A T 

2 rria 



TV? 



X B 

2tt'^ 771^ ' 



(23) 



Where A is the t'Hooft coupling, T is the temperature and B and are the magnetic field and 
the bare mass of the dual gauge theory. The phase diagram of the theory in these coordinates 
is presented in figure |4} 



The fitting curve for large B/77? corresponds to a 



behaviour. This refiects 



the existence of a critical ratio of the magntic field and the temperature at witch the phase 
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confined phase 




Figure 3: A phase diagram of tlie tlieory. Tlie blue curve separates tlie confined and deconfined phases 
of the theory. There is a first order phase transition across the critical curve. For 77 < r/cr ~ 4.25 
at zero bare mass the theory has a non-broken global S0(3) symmetry, while for rj > rjcr the global 
symmetry is spontaneously broken and the theory has a mass-gap. 



transition disappears rjcr- Indeed expanding m{ri) near rjcr and using that m(?7c 
Therefore to leading order we have: 



1 

fh 



V 



we obtain: 
(24) 

(25) 



which given the relations from equation (23) is equivalent to the observed: 



2 m„ 




const; const 



(26) 



behaviour for large B/m'^. On the other side the fitting curve for small B /m'^ in figure |4] 



\ B \2 



behaviour. This simply reflects the fact that /ds is a function 



corresponds to a ~ 
of even powers of r] as evident from the definition in equation (|21j) 



It is interesting to compare our phase diagram to similar phase diagrams for 1 + 2 dimensional 
field theories in both external magnetic field and finite temperature. It is somewhat intriguing 
that the phase diagram of the 1 + 2 dimensional Gross-Navie model studied in refs p!7l [18] 
has the same qualitative structure. It would be interesting to use alternative non-perturbative 
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techniques to study the phase diagram of the defect field theory holographically dual to our 
set up and compare with the result obtained via the AdS / CFT correspondence. We leave such 
studies for future investigations. 



3.3 The order parameter 

In this subsection we focus on the zero bare mass case when the confinement /deconfinement 
phase transition is also a spontaneous symmetry breaking phase transition. The quantity that 
we are interested in is the fundamental condensate of the theory which serves as an order 
parameter of the transition. In particular we are interested in the temperature dependence of 
the condensate at zero bare mass. 

Let us consider again the zero temperature case with no external field. As mentioned above 
at zero bare mass the theory is conformal and has a global 5*0 (3) x 5*0(3) symmetry. If we 
introduce external magnetic field [31] the conformal symmetry is broken and the theory has 
meson-like excitations. Furthermore the theory develops a negative fundamental condensate 
that breaks the S0{3) part of the global symmetry down to a U{1) symmetry. There are also 
zero modes identified with the goldstone bosons of the spontaneously broken global symmetry. 
On the other side if we turn on a finite temperature the adjoined degrees of freedom are in a 
deconfined phase. Our holographic set up now describes pure Yang-Mills plasma with some 
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dilute fundamental matter. At low temperatures the fundamental matter is still in the confined 
phase and is better described as a gas of mesons. However, at sufficiently high temperature the 
alternative deconfined phase becomes the stable one and the theory undergoes a first order phase 
transition. The mesons are "melted" and the dynamically generated mass of the fundamental 
hypermultiplet fields as well as their condensate vanish. In the holographic set up this is 
refiected by the fact that the / = embedding of the D5-brane becomes the stable one. This 
suggests that the global symmetry of the theory is restored in the deconfined phase and we can 
reffer to the transition "chiral" transition. 

The physical picture described above can be visuahzed if one generates a plot of the temperature 
dependence of the fundamental condensate (for a fixed magnetic field). More precisely we 
generate a plot of the ratio of the condensate at finite temperature and the condensate at zero 
temperature {il!ilj)T/{4''ip)o versus the ratio of temperature and the square root of the magnetic 



field T jyB. To this end we use that according to equations ([9j) and (15) we have the relations 



(V'V')t 1 . 
oc — c; 



B 



T 



(27) 



Note also that the ratio B is obtained in the rj 

plot is presented in figure [5] 



oo limit in equation (27). The resulting 
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Figure 5: A plot of {'iIji{j)t/ {'ip'4^)o versus T/VB. The vertical dashe line corresponds to the critical 
temperature at which the "chiral" transition takes place. The purple colour corresponds to either stable 
or metastable phases while the red colour represents unstable phases with "tachyonic" excitations. The 



fitting curve at small temperatures represents the 1 — const T'^ dependence from equation (28). 
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The upper purple curve corresponds to the confined symmetry breaking phase, the colour 
represents that this phase is "tachyon" free. (As we shall show in Section 4). The states on 
the left-hand side of the vertical dashed line are thermodynamically favoured and correspond 
to a stable phase of the theory. The states on the right-hand side of the critical curve are 
metastable states. The red curve with negative slope corresponds to unstable phase of the 
theory which poses a tachyon in the spectrum. It is also unstable from thermodynamic point 
of view. The horizontal red curve corresponds to a symmetric deconfined phase of the theory 
which is unstable for temperatures below the critical value (represented by the vertical dashed 
line). The purple horizontal line corresponds to a symmetric deconfined phase of the theory 
which is "tachyon" free. The states on the right-hand side of the vertical dashed line are 
the thermodynamically favoured ones and correspond to a stable phase of the theory. The 
states of the horizontal purple line that are on left-hand side of the critical line correspond to 
a metastable phase of theory and can be reached by supercooling. The fitting curve at low 
temperatures represents the asymptotic relation: 

||l = l-coust(— j^; const. 0.929; (28) 

This relation can be obtained by expanding the equation of motion of the D5-brane embed- 
dings in powers of l/r/ and solving perturbatively for the leading corrections of the D5-brane 
embedding. The perturbative solution is correction near the zero temperature solution ob- 



tained in [31]. It would be nice to derive equation (28) from field theory side. One can also 
study this relation in the context of chiral perturbation theory. We leave such studies for future 
investigations. 

3.4 The entropy denisty 

In order to calculate the entropy density of the theory and verify its positivity we can use 



equation (22) for the free energy density. The entropy is given given by: 



= _Nf^An'R%' ( 3/d5 + mc- = Nf^An^R^'^Sim, r/^); , 

where we have used that c = —9^/05- We have also defined the dimensionless quantity S{m, rf) 
which is convenient to study numerically. 
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It is instructive to calculate the entropy of the /(p) = embedding first. That is the entropy 
at zero bare mass in the deconfined phase with non-broken global symmetry characterized by 
vanishing condensate. It turns out that in this case the entropy can be obtained in closed form. 
The result is: 

5(0, = yrr^ . (30) 

In field theory units the entropy is: 



S{mg = 0, T, B) = ^/Ij2NfN,\jT^ + — . (31) 

Note that at large temperature as T ^ oo (or equivalently weak magnetic field) the entropy 
grows as oc as one would expect for a 1 + 2 dimensional conformal field theory. However 
the limit of zero temperature (at fixed magnetic field) doesn't lead to vanishing entropy. This 
apparently strange behaviour is rectified if one takes into account the phase diagram of the 
theory. Indeed, at fixed magnetic field the limit of zero temperature suggests that ?7 — oo 
which is clearly above the critical value of rjcr and the stable phase of the theory is the confined 
phase with spontaneously broken global symmetry. Therefore for very small temperature the 
entropy at zero bare mass should be evaluated at the state characterized by both vanishing 
bare mass and (the smallest) negative condensate. The intermediate picture has a first order 
phase transition aX t] = with a finite jump of the entropy corresponding to the latent heat 
of the confinement /deconfinement phase transition. At zero bare mass this transition is also 
a "chiral" phase transition. In fact at small (and zero bare mass) the jump of the entropy is 
enhanced because the dynamical mass generation associated to the "chiral" phase transition 
freezes some part of the light degrees of freedom (in the symmetry breaking phase). The latent 
heat is further increased (relative to the zero magnetic field case) because of the increased 
diamagnetic response of the theory in the deconfined phase (look at subsection 3.4). 

To generate a plot of the entropy S versus the temperature T at fixed magnetic field B and 
bare mass m^, we use the following relations: 

1 fXi^^V T 71 S 1~ 

S(m,r/); (32) 



^ V 2 ; v^' NfN.B r] 

Therefore to evaluate the entropy at zero bare mass we need to generate a plot of ^5'(0, t]) versus 

-p. The resulting plot is presented in figure 6 The red curve in the figure corresponds to the 
v'' I I 

deconfined phase of the theory. As on can see near the phase transition the entropy is a multi- 
valued function of the temperature. The upper smooth red curve corresponds to a phase with 
vanishing condensate and non-broken 5*0(3) symmetry. It is generated using equation (30) for 
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Figure 6: A plot of the entropy for zero bare mass niq = 0. The upper red curve is generated using 
equation (30) for S{0,r]). The red curve with negative slope is generated numerically and corresponds 
to an unstable deconfined phase with broken global symmetry. The blue curve corresponds to a 
confined phase with broken global symmetry. The fitting curve at large temperature corresponds 
to dependence as expected for 1 + 2 dimensional conformal theory. The fitting curve for small 
temperature corresponds to behaviour. There is a first order "chiral" phase transition at ry = r]cr 
(the vertical dashed line) characterized by a finite positive jump of the entropy denisty. Interestingly 
the ratio of the entropy immediately before and after the phase transition is equal to 27r. 



S'(0, 77) . The branch with negative slope corresponds to a unstable phase characterized by non- 
zero condensate and is generated numerically scanning through the corresponding Black hole 
D5-brane embeddings. The lowest lying blue curve corresponds to a confined phase in which 
the global symmetry is spontaneously broken. It is generated numerically and corresponds to 
Minkowski D5-brane embeddings. At rj = rjcr there is a first order phase transition (represented 
by the vertical dashed line in the figure). One can see that there is a positive jump of the entropy 
as going from the confined to the deconfined phase. The fitting curve at high temperature 
corresponds to behaviour, as expected from a 1 + 2 dimensional conformal field theory |13j 
(the effect of the external magnetic field is suppressed in this limit). At low temperature one 
can show that the entropy vanish as T^. The coefficient of proportionality can be calculated 
numerically. 

For generic bare masses the entropy can be studied in two different regimes: the regime of fixed 
ratio of magnetic field and temperature (fixed r]) and the regime of fixed ratio of bare mass 
and external magnetic field (fixed nig/^/B). The former regime is technically more convenient 



19 



to study (the quantity Ids depends explicitly on 77). This is also the regime studied in refs. 
j25l [21]. However in this regime the slope of the entropy 5* versus temperature T plot is not 
proportional to the specific heat at fixed magnetic field cb, because the magnetic field varies 
with temperature (in order to keep 77 oc B/T"^ fixed). This is why we study the entropy at 
fixed ratio of the bare mass and the magnetic field. This is also a natural generalization of the 
study of the zero bare mass case and as we are going to see the qualitative behaviour remains 



unchanged. To perform the study we use the dimensionless variables defined in equation (32). 
The parameter that we keep fixed is: 

The resulting plots for m = 0.5, 2 are presented in figure [7| Similarly to the zero bare mass 





Figure 7: Plots of entropy versus temperature for m = 0.5,2. The qualitative behaviour is similar to 
that from figure |6j Major difference is that the deconfined phase with non-broken global symmetry 
does not exist for temperatures smaller than the temperature at which the negative slope deconfined 
phase appears. 

case at high temperatures the entropy grows as which is the dependence expected for 1 + 2 
dimensional conformal field theory at finite temperature. For intermediate temperatures there 
is a finite jump of the entropy associated to the confinement/deconfinement phase trasition of 
the fundamental matter. At low temperatures the entropy goes to zero as T^. The coefficient of 
proportionality depends on the ratio m^/ ^/B and can be calculated numerically. Note that near 
the phase transition the entropy is a multi-valued function of the temperature. The multiple 
phases regime of the theory exists for a finite interval of temperatures and at the boundaries 
of this interval the specific heat cb diverges. A more detailed study of the entropy function in 
the unstable phase would reveal that the negative slope curve has more complicated structure. 
Namely, there is a double logarithmic structure and the red and blue curves join in the centre 
of that spiral [7j. We shall not study further this property of the theory. 
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To determine the critical value of 1/ ^ we have used the phase diagram from figure |4| Al- 
ternatively we could have used the equal areas law |25], because the area below the curves in 
figure [7] is proportional to the free energy of the system. 



3.5 The magnetization 

In the statistical ensemble that we consider, that is fixed temperature T and fixed magnetic 
field B, there is another natural quantity of interest. This is the magnetization of the theory 
given by: 

M^-m ^-S.^.-N,t,,^m f - (34) 



\dB J rp Qs \ drj j W' 



m 



Here we have defined the dimensionless quantity M. In filed theory units the last expression 



in equation (34) can be rewritten as: 



M = -^TM . (35) 

Thus the quantity M is convenient for studies at fixed non-zero temperature. In order to 
incorporate the zero temperature case it is convenient to define: 



where we have used the first equation from (32). The dimensionless quantity M is given by: 



oo 



pmin 

Note that we have introduced a new set of dimensionless variables: 



p = p/rVH; i = 1/rVH; (38) 

convenient to study the theory at a fixed magnetic field. It is instructive to study first the 
magnetization at zero temperature. 
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3.5.1 Magnetization at zero temperature 



The zero temperature case corresponds to the hmit 1] —>■ oo in equation (37). The expression 
for the magnetization M is given by: 

oo / ^ — 

M„ = - / d /^f±g . (39) 

Pmin 

Here l{p) is a solution to the equation of motion for the D5 brane probing a pure AdS^ x 
space. We shall be interested only on the stable branch of the theory. For large bare masses 
m ^ 1 one can use the results from ref. [31] for the asymptotic behaviour of the condensate: 

■ ' (40) 



to show that the asymptotic behaviour of the magnetization should be: 



M„ = -|-. (41) 



after using equations (33) and (36) we obtain: 



NfN B 

Therefore at weak magnetic field (large bare mass) the meson gas has diamagnetic properties 
with negative susceptibility: 

dMo NfN, 1 

Xo = -WW- = i. + . . . . 43 

oB An uig 

For large magnetic field the magnetization can be obtained numerically. The resulting plot is 



presented in figure [8j As one can see equation (41) is a good approximation for large bare 
mass (weak magnetic field) while at zero bare mass the magnetization has a non-vanishing 
negative value. Furthermore the finite positive slope at small rh suggests finite negative sus- 
ceptibility. This confirms that at zero temperature the meson gas is in a diamagnetic phase. 
One may expect that at finite temperature the diamagnetic properties of the meson gas would 
not depend strongly on the temperature. However at high enough temperature a confine- 
ment/deconfinement phase transition takes place and the theory is in a conductive phase. It 
would be interesting to study the effect that this transition has on the magnetic properties of 
the matter. 
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Figure 8: A plot of Mq oc Mq/VB versus m oc niq/yB. The dashed curve corresponds to equation 



(41) one can see the good agreement for large bare masses. At zero bare mass the magnetization has 



a finite negative value. The finite positive slope at small m suggests finite negative susceptibility. 
3.5.2 Magnetization at finite temperature 



It is instructive to study first the magnetization at zero bare mass = 0. At zero bare mass 
the theory can be in either confined or deconfined phase depending on the ratio of the magnetic 
field and the temperature squared. This is controlled by the parameter rj. For small values 
of r] the theory is in the deconfined phase and has a non-broken SU(3) global symmetry. The 
corresponding D5 brane embedding is the / = one. For large enough temperatures (weak 
magnetic field) M behaves as: 

M = -r] + 0{rf) . (44) 
This suggests the following expression for the magnetization M in field theory units: 



M 



1 



O I ^ 1 



T2 J 



(45) 



Clearly the magnetization is negative and the leading contribution to the magnetic susceptibility 



is: 



X 



2^/'^NfNc 1 
AV2 



T 



(46) 



Therefore the deconfined phase of the theory is a diamagnetic phase. Furthermore the dia- 
magnetic response depends strongly on the temperature and goes to zero as the temperature 



approaches infinity. Curiously equation (46) is similar to the Curie's law % oc 1/T for a para- 
magnetic but with the wrong sign. This is also the most simple expression for the magnetic 
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susceptibility in 1 + 2 dimensions based on dimensional analysis only. Such a behaviour is to 
be expected since at high temperatures conformality is restored. 

For strong magnetic field the magnetization can be obtained numerically. It is convenient to 



use the parameter M oc MjyB defined in equation(36). A plot of M versus l/y^ is pre- 



sented in figure [9) Note that the axes are labeled in field theory units. The blue curve for low 
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Figure 9: A plot of M versus l/^/rj. The axes are labeled in field theory units. The blue curve 
corresponds to the confined global symmetry breaking phase of the theory. The dashed horizontal line 
corresponds to Mo|rn=o- The dashed vertical line denotes the first order confinement/deconfinement 
phase transition. The fitting curve at high T corresponds to equation (|45|). 



temperature corresponds to the confined global symmetry breaking phase of the theory. The 
dashed horizontal line corresponds to Mo|m=o- As one can see the magnetization is negative 
and the diamagnetic response of the theory varies slightly with the temperature. At sufficiently 
high temperature (the vertical line in the figure) the theory undergoes a first order confine- 
ment/deconfinement phase transition, which is also a "chiral" phase transition. As one can see 
there is a significant jump of the magnetization suggesting much stronger diamagnetic response 
of the theory in the deconfined phase. This is expected because the confined phase is also a 
conductive phase. Upon further increase of the temperature the magnetization approaches the 



behaviour described by equation (45) (the fitting curve for high T). One can see that in the 
deconfined phase the magnetization depends strongly on the temperature. One explanation is 
that the conductivity drops as the temperature is increased. 

Another interesting feature is that for the minimal value of 1/ y/fj for which the deconfined 
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exists there is a kink and slope has a finite jump. This is in contrast to the point corresponding 
to the maximum value of 1/ for which the confined phase exists. The slope diverges at this 
point and the blue curve bends smoothly. It seems that this is specific for the rriq = case. 

The case of generic values of rh can also be studied numerically. Plots of M versus 1 / ^/f] for 
TTi = 0.5, 2 are presented in figure [T0| As one may expect the magnetization is negative and for 
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Figure 10: Plots of M versus for rh = 0.5,2. The blue curves correspond to the confined phase 

and the horizontal dashed lines correspond to Mo|m=o.5 and Mo|rft=2 correspondingly. The vertical 
lines correspond to r/cr for which the confinemet /deconfinement phase transition takes place. The 
red curves describe the magnetization in the deconfined phase of the theory and the fitting curves 



correspond to equation ( 45 1 . 



temperatures below the critical one the theory is in a confined phase (the blue curves in the 
figure). One can see that the magnetization varies slowly with temperature and approaches very 
fast the values for vanishing temperature (Mo|m=o.5 and Mo|rn=2) represented by the horizontal 
dashed lines. At the confinement /deconfinement transition there is a jump of the magnetization. 
It is interesting that the relative jump of the magnetization seems to decrease as rh increase 
(m = 0, 0.5, 2). For high temperature the theory is in a deconfined phase and the magnetization 
approaches the behviour described by equation (45). Unlike the m = case there is no kink at 



the point where the deconfined phase cease existence. 



4 Meson Spectrum 

In this section we study the light meson spectrum of the theory. The holographic prescription 
for obtaining the spectrum of meson-like excitations is to study semi-classical fiuctuations of the 
probe fiavour branes. Let us consider a small excitation of the probe brane propagating along 
the radial direction of the geometry. Since the geometry is asymptotically AdS space-time the 
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excitations would bounce at infinity without dissipation. On the other side the geometry has a 
horizon in the bulk and there are two possibilities. If the D5-brane has a shrinking S*^ and close 
above the horizon (Minkowski type of embeddings) the excitations would bounce and as a result 
the D5-brane would have normal modes corresponding to standing waves. Another possibility 
is that the D5-brane reaches the horizon before the S"^ shrinks (Black hole embeddings) in 
this case the excitation fall into the horizon and dissipate. The corresponding semiclassical 
excitations are quasi-normal modes satisfying ingoing boundary conditions at the horizon. 

The interpretation on gauge theory side of the correspondence is that the normal modes of 
Minkowski embeddings correspond to the meson spectrum of the confined phase, while Black 
hole embeddings correspond to the deconfined phase of theory and their quasi-normal modes 
describe mesons melting in the Yang-Mills plasma [12]. The real components of the quasi- 
normal modes correspond to the energy of the melting meson and the imaginary component 
correspond to the relaxation time (if it is negative). 

A study of the semiclassical excitations of the theory is of a particular interest when we have 
multiple phases because it verifies the thermodynamic analysis of the stability of different 
phases (existence or non-existence of tachyonic modes). It is also important when we have a 
spontaneous breaking of a continuous symmetry because of the existence of Goldstone modes 
and related phenomena. 

The meson spectrum of the holographic gauge theory corresponding to the D3/D5 set up has 
been extensively studied in the literature. In refs. [39l HO] the spectrum of the supersymmetric 
zero temperature case has been investigated. In ref. an external magnetic field has been 
introduced leading to a non-supersymmetric field theory exhibiting spontaneous symmetry 
breaking. The corresponding meson spectrum has been investigated with a focus on the pseudo- 
Golstone modes. A non-relativistic dispersion relation of the goldstone modes has been revealed 
and various phenomenological relations verified. However analysis of the stability based on the 
massive meson modes has not been performed. In this paper we shall perform such studies in 
the more general case when a finite temperature is turned on. 

In order to study the light meson spectrum of the theory we look for the quadratic fluctuations 
of the D5-brane embedding along the transverse directions parametrized by /, tp, cf). To this end 
we expand: 

1 = 1 + 27ra'6l; = 27ia'6ij; = 27ra'6(f) , (47) 

in the action ^ and leave only terms of order (27ra')^. Note that fluctuations of the U{1) gauge 
field Faf3 of the D5-brane will also contribute to the expansion. There is also an additional 
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contribution from the Wess-Zumino term of the D5-brane's action: 



5wz = NffX5 [Yl^CpAe^]; T^B + 2na'F . (48) 
Me P 

For the anzatz that we are considering, the relevant term is: 

Swz^Nfi^^ J BAP[C4 , (49) 

Me 

where -PfGi] is the pull-back of the magnetic dual, C4, to the background C4 R-R form. For 
the particular parameterization of 5*5 considered here, it is given by: 

1 Ao'^P 

C4 = - . / ^^. R^sinipUdp - pdl) A (1^2 A # . (50) 
Qs [p^ + i^r 

After some long but straightforward calculations we get the following action for the quadratic 
fluctuations along I: 



5f , (51) 



'(2) _ ^^^ ^a 712 _ Pi 7127'^ 
1- 



4 

and along and -0: 

°^ ^v^5"^(G'^^a„<5V'9/35V' + G.^^a^^c/.a^j^^) , (52) 
C^^^ oc {cos a) PHSipdoScj) . 
Here E'^/j is the pull-back of the generalized metric on the classical D5-brane embedding: 

Eaf3 = d^X^dfiX^G^, + B^,) , (53) 

and we have defined 8°"^ and 3°"^ as the symmetric and anti-symmetric elements of the inverse 
generalized metric E""^: 

The functions P and 5f(p) = sJ—Ej cos a are given by: 

4d4 2/2 
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4.1 Fluctuations along 1. 



We first focus our attention to tfie fluctuations along I. Note that there is no symmetry 
associated to translations along / and we do not expect the appearance of Goldstone bosons. 
However the absence or presence of tachyons in the spectrum of 61 can give is information about 
the stability of the phase that we are studying. 



The equation of motion derived from the quadratic action (52) is the following: 



;i + /'2)s 



+ 



;i + /'2)(4r4 + 64) 



B5l + ^^5l- 



dfgip) 



d 

dp 



l'di9{p)\ 



where 



□ 



-(92 



+ 



1 + 



51 = 0(57) 



(58) 



and A (2) is the spherical laplacian. Note that the modified operator of d'Alembert reflects the 
breaking of the S0{1,2) Lorentz symmetry down to a S0{2) symmetry. 

Let us focus on the meson spectrum of the various phases at zero bare mass. In particular let 
us verify that the symmetry breaking phase corresponding to a confined phase of the theory 
(the upper purple curve in figure [s] ) is stable or metastable and has no tachyon modes. To this 
end we consider a plane- wave ansatz: 



Sl = ap)e 



—iujt+ik.x 



(59) 



Since we are interested only on the energy of the excitations we will further restrict our anzats 
to A; = 0. the equation of motion can be written as: 



'i+i 



i2\2 



4^(^)(4f^+l^ 



;i + /'2)(4f 



?:4 



dtm - 



d ( l'dfg{p) 
dp\l + /'2 



C = o. 



(60) 



Where we have rewritten the equation in dimensionless variables and defined uj = R^uj/h. 



4.1.1 Normal modes. 



Next we consider the Minkowski embeddings that asymptote to zero separation m = at 
infinity and solve numerically the equation of motion (60) imposing the following boundary 
conditions at p = e = 10~^: 

C(e) = l; C'(e) = 0; (61) 



the spectrum is obtained by requiring regularity of the solution at infinity. 
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To study the ground state of the meson spectrum in the symmetry broken phase and its 
temperature dependence we generate a plot of cu/cjo versus T/\fB using that: 



u 1 

OC ^—OJ 



and that the zero temperature value of the frequency Uq corresponds to the r] ^ oo limit. A 
plot of of uj/ujo versus l/y/rj is presented in figure 
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Figure 11: A plot of the spectrum of the ground state of the symmetry breaking phase. The spectrum 
remains real even for temperature above the critical temperature, represented by the vertical dashed 
line. The spectrum drops to zero at the state corresponding to divergent heat capacity. 



As one can see for low temperature the ground state is tachyon free and the confined symmetry 
breaking phase is the stable one. Furthermore even for temperatures above the critical one 
(the vertical dashed line in the figure) the symmetry breaking phase remains metastable. If we 
further increase the temperature the energy of the ground state decreases sharply. Remarkably 
at the maximal temperature for witch confined phase exists the energy of the ground state 
drops to zero. Note that at this point the heat capacity also diverges as the plot of the entropy 
from figure [6] infers. Note that near the maximal temperature there is an alternative confined 
phase represented by a red curve in figure [5] and having negative slope in figure |6} A study of 
the meson spectrum shows that this phase has tachyonic ground state and is unstable. The 



spectrum is represented by the blue curve in figure 12 
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4.1.2 Quasi-normal modes 



To study the deconfined phase we need to study the qusi-normal excitations of the probe D5- 
branes. This suggests imposing ingoing boundary conditions at the horizon of the geometry 
this is why it is more convenient to use spherical coordinates in the transverse subspace. 
The equations are further simplified if we introduce a new radial coordinate: 

2 4r^ + 6^ 



u 



In the new coordinates the metric of the AdS^ x black hole is given by: 



ds^ /a 



u 



u 



dt + — rrff^ + 



rdu^ + R\de'' + cosO'^d^i + sin^2f]2^ 



(63) 



(64) 



^2 - ' m4_^4- 

note that in these coordinates fluctuations along I corresponds to fluctuations along 9. It is also 
a straightforward exercise to rewrite the equation of motion for 51 ( [57| ) in terms of 6{u) and 
59{u). It is convenient to introduce dimensionless coordinates u = u/b. The resulting equation 
of motion is given by: 



g{u){u^ - 1)56' 



dgg{u) - du 



+ 



9W 



056 + 



9[u) 



A 



(2) 



deg{u) I 56 = 



K+^6'^ cos ^2 



56 



(65) 



1 + ^6'^ 



where 



9[U) 
□ = 



cos 6^^ 



1 + ^ 



2\ 1/2 



r,4 



1 + 



u 



1/2 



(66) 
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62 



-dl ^&i+d, 



1-^ 
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To study the spectrum of the fluctuations we consider an anzats: 

56 = e-'"*C(M) 



(67) 



Note that in general is a complex number. The real part is naturally interpreted as the energy 
of the excitation, while the imaginary part of u is proportional to the relaxation time of the 
excitation, provided that it is negative (Im(u;) < 0). Clearly a positive imaginary part would 
lead to exponentially growing (with time) excitations. The existence of such modes signals 
instability of the phase under consideration. 
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To quantize the spectrum we need to impose ingoing boundary conditions at the horizon. To 
this end we focus on the equation of motion for ({u) near u = 1. To leading order we have: 

C'iu) + ^C{u) + TT^T^^A^) = , (68) 
u — 1 lb[u — 1)^ 



where we have defined u = ^uj. The most general solution of equation (68) is a linear 
combination of {u — i)='=*'^/4. However ingoing boundary conditions correspond to the negative 
sign solution {u — i)^*'^/^ (^for the choice of signs considered in (67)). To impose the condition 
we define: 

C,{u) = {u-\r^l'S{u)- (69) 

and solve numerically the resulting equation of motion for S{%l) imposing Dirichlet boundary 
condition at the horizon (5'(1) = 1;). The quasi-normal modes are obtained by selecting regular 
solutions at infinity. More precisely we require that \{u — l)~*'^/^S'('u)| ~ l/tt^ as we send 
u ^ oo. Next we consider the classical black hole embeddings which have a zero separation 
at infinity m = 0. There are two classes of solutions. Solutions which develop separation in 
the bulk of the geometry and hence describe deconfined symmetry breaking phase and the 
6' = solution which has restored symmetry. It turns out that that the spectrum of quasi- 
normla modes for both type of solutions contain a mode which is purely imaginary and thus 
particularly convenient to study numerically. A plot of the spectrum is presented in figure [12} 
The smooth red curve corresponds to the ^ = solution and thus describe a deconfined phase 
with restored global symmetry. As one can see for large temperatures the imaginary part of 
the quasi-normal mode is negative and thus the phase is stable. As we lower the temperature 
we reach the critical temperature represented by the vertical dashed line. At this point a first 
order phase transition takes place and the thermodynamically favoured phase is the confined 



symmetry breaking phase. However as one can see from figure 12 the symmetric phase remains 
metastable even for temperatures below the critical one and thus describe states that can be 
reached by supercooling. If we continue to lower the temperature the imaginary part of the 
quasi-normal mode becomes positive and the deconfined phase is unstable. Remarkably this 
happens at the point where the symmetry breaking black hole solutions and the 6* = solutions 
meet. In the limit T ^ the spectrum approaches the spectrum for pure AdS'^ x space 
represented by the horizontal dashed line. 

The symmetry breaking deconfined phase is represented by the short red curve with positive 



slope in figure 12 (and negative slope in figure |6]). As one can see it has quasi-normal excitations 
with positive imaginary components and is thus unstable. This fits nicely with the fact that it 
has a negative heat capacity (negative slope in figure [6]). 
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Figure 12: A plot of the imaginary part of the of the quasi-normal modes (the red curves) and of 
the tachyonic sector of the normal modes (the blue curve). The vertical dashed line corresponds to 
the critical temperature at which the confinement /deconfinement phase transition takes place. The 
horizontal dashed line corresponds to the highest quasi-normal mode at zero temperature. 



Finally the blue curve in figure 12 corresponds to Minkowski solutions which have tachyonic 
meson spectrum and is obtained by solving equation (60). This phase is represented by a blue 



curve in the entropy versus temperature plot (figure [6]) and has a negative slope. Clearly it 
is unstable as both the presence of tachyons and negative heat capacity suggest. Note that 



the kink in figure 12, where the red and blue curves join, corresponds to embeddings near the 
critical embedding which separates black hole embeddings and Minkowski embeddings. A more 
detailed analysis of the theory in this regime would reveal a discrete self-similar behaviour [13] . 
We shell not elaborate further on this property of the set up. 



4.1.3 The zero magnetic field case 

For completeness and to verify the validity of our numerical analysis it is instructive to study 
in details the spectrum of quasi-normal modes at finite temperature and zero magnetic field. 
At zero bare mass we can employ the technique used in ref. [12] and solve the corresponding 
Heun equation. 



To begin with let us write down the equation of motion ( 66 ) for the fluctuations of the 9 = 
embedding: 

+ - :r^) mu) + ( ^ + SOiu) = (70) 



u{u^ + ri'^)J V(m4-1)2 
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Next we focus on the 77 = case and consider the substitution x = 1 — l/u"^. The equation of 
motion is written as: 

•■ l + 3(l-x)2 



2x{l - x){2 - X 



-69 + 



+ 



4x2(1 - 2;)(2 -a;)2 4x(l - x)2(2 - ,t 



69 = 



(71) 



Equation (71) has four regular singularities x = 0,1,2, 00 with exponents {—iw/A,+iw/A}, 
{1/2,1}, {—u/4:,+Lj/4:}, {0,0}. Upon the change of variables: 



69{t,x) 



-iLjt„—i 



X *4 (1 — x){x — 2) 4y(x) 



equation (|71|) takes the standard form of a Heun equation 

a/3-g 



y{x) + ( - + 3^ + ) y[x) + 



X X — 1 X — 2 



x{x — l)(x — 2) 



y{x) = , 



(72) 



(73) 



with parameters: 



^ = l-iw/2; 6 = 3/2; e = l-Cu/2; Q = 3/2 - (l/A + i)Cu - {l/A - i/8)Cu^; (74) 
a = 3/2 - (1/4 + i/A)u; P = l- (1/4 + i/A)^; e = a + p~-f-6 + l; . 



Next we look for solutions of equation (73) satisfying y{0) = y{l) = 1. This selects solutions of 



equation (70) obeying ingoing boundary condition at the horizon {x = 0) and an appropriate 



behaviour at infinity (x = 1). The method that we use is the one employed in ref. [I2] for the 
D3/D7 set up. We consider a Frobenius series near x = satisfying the recursion relation: 



a„+2 + An{uj)an+i + B{uj)an = , 



where An{uj) and Bn{uj) are given by: 



(n + l)(2(5 + /t + 3(n + 7)) + g 



2(n + 2)(n + l + 7) 
n + a){n + (3) 



^"^""^ 2(n + 2)(n + l + 7) 



(75) 

(76) 
(77) 



and = Q/2'j. Next we define a continued fraction [38], [T2] : 



BJCj) 



AJuj) +r„+i 



(78) 



the convergency condition is given by: 



ro = Q/27 • 



(79) 
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To obtain approximate expression for tq we cut the recursive relation at some sufficiently large 
n {n = 150 in our case) and use the asymptotic expression for r„: 

_ 1 2 + Co 



+ 



(80) 



Next we solve the resulting algebraic equation for uj. The spectrum for the first 11 lowest 
quasi-normal modes is given in table [Tj Qualitatively the spectrum of quasi-normal modes has 





k 


Re{Co) 


Im(ti;) 





±1.6906 


-1.3268 


1 


±3.6805 


-3.3269 


2 


±5.6780 


-5.3261 


3 


±7.6769 


-7.3258 


4 


±9.6764 


-9.3256 


5 


±11.676 


-11.326 


6 


±13.676- 


13.325 


7 


±15.676 


-15.325 


8 


±17.676 


-17.325 


9 


±19.675 


-19.325 


10 


±21.675 


-21.325 



'3 



■10 



■15 



-20 



-20 



-10 





Reicj) 



10 



20 



Table 1: Quasi-normal modes for the first 11 excited states. 

the same structure as the one for the D3/D7 system studied in ref. [T2j. As one can see from 
the plot in table [T] for large k the quasi-normal modes satisfy the approximate relation: 



Uk = ±0.68 - .33i ± {2k + 1)(±1 



(81) 



If we apply the numerical techniques used to generate the spectrum in figure 12 to obtain the 
spectrum presented in table [l] we find that the results agree with relative error of 0.01% (for 
the lowest quasi- normal modes). A more complete and detailed study of the spectrum of quas- 
normal modes for non-zero magnetic field is beyond the scope of this paper. However analysis 



of equation (70) shows that as the magnetic field increases the lowest lying modes given in 
table [l] shift radially away from the origin. Furthermore a new set of quasi-normal modes with 
vanishing real part emerges. The lowest lying of these modes is the one analyzed in the previous 



subsection. As one can see from figure 12 for sufficiently large magnetic fields this mode becomes 
tachyonic. If we keep on increasing the magnetic field more and more tachyonic modes emerge. 
In the strict r] —>■ oo limit (corresponding to vanishing temperature) the spectrum contains an 
infinite tower of tachyonic modes. The qualitative behaviour is similar to the one described in 
where the zero temperature spectrum for the D3/D7 set up was explored. 
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4.2 Fluctuations along and tjj 



In this subsection we study fluctuations along and ijj. Since translations along (f) and ^p 
correspond to the generators of the spontaneously broken SU (2) symmetry for small bare 
masses we expect to detect pseudo-Goldstone modes. This is why we focus on the spectrum 
in the confined phase. The spectrum at zero temperature has been extensively studied in ref. 
|31j . We expect that the qualitative behaviour at finite temperature will remain the same in 
the confined phase. 



The equations of motion derived from the quadratic action (52) are the following: 



4;2 



d 



dM I + 

+ 



9{pr 



A(2)(5^ - PHdo 







(82) 



066 + 



P' 



where □ is give in equation (58). Next we consider a plane- wave ansatz: 



r/i(p); 5i) 



(83) 



now using the anzatz (83) we get: 

'9ipr 



1 +/'2 

9{pr 

1 + /'2 ■ 



4/2 



^2 + 



^9{p)R'l 
4r2 + 64 

^9{p)R'l^ 
4^4 ^ 54 





^2 




P 


1 - 




1 + 




(4r4+fe4)2 


(4r4+64)2 
P 


1 - 


1664r4 


1 + 


16R4H2r4 


(4r4+f)4)2 


(4r-4+64)2 



r/i - iujPHri2 = , (84) 



r]2 + iuPHrji = 



The equations of motion in (84) can be decoupled by the definition ri± = rji ± ir]2. The result 
is: 



9{pr 

1 +/'2 

9{pr 

1 +/'2 



^9{p)R'l^ 
4^4 + 54 

^9{p)R'l' 
4r4 + 64 



V- 





u;2 




P 


1 - 


16b4r4 


1 + 




(4r4+fe4)2 

a;2 


(4r-4+64)2 
P 


1 - 


16b4r4 


1 + 




(4r4+fe4)2 


(4r4+b4)2 



?7+ — uPHt]^ 







^5) 



r]^ + luPHt]^ = 



We can now solve numerically the equations of motion by imposing boundary conditions: 



(86) 



For sufficiently small e ~ 10^6. The spectrum is quantized by requiring that r^^ ~ 1/p for large 
p. Since we are interested in the confined phase we consider 77 = 4.5 > rjcr- The resulting spec- 
trum is presented in figure 13, where we have introduced dimensionless parameters m = m/b 
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Figure 13: A plot of tD versus m. One can observe a Zeeman splitting of the energy levels at large 
bare masses m. And the existence of a pseudo-Goldstone mode at small bare masses. 



and w = B?uj j b. As one can see for large bare masses there is a Zeeman splitting of the energy 
levels, and the spectrum approximates to the spectrum for zero temperature and magnetic field 
studied in refs. [311 HO], where the authors obtained the following relation: 

cc'„ = |^V(ri+l/2)(n + 3/2) , (87) 

between the eigenvalue of the n^^ excited state Un and the bare mass m. The black fitting lines 



in figure 13 correspond to equation (87). It is also evident that at zero bare mass the ground 
state has zero frequency. This is the Goldstone mode of the broken SU (2) symmetry. Note 
that we have broken two generators while we observe only a single Goldstone mode. This is the 
same behaviour as in the zero temperature case studied in ref . [21] . The apparent contradiction 
is clarified by the observation that the S0(l,2) Lorentz symmetry is broken down to S0(2) 
rotational symmetry (by both the external magnetic field and the finite temperature). This 
opens the possibility of having two types of Goldstone modes: type I and type II satisfying odd 
and even dispersion relations correspondingly. In this case there is a modified counting rule 
(ref. [H]) which states that the number of GBs of type I plus twice the number of GBs of type 
II is greater than or equal to the number of broken generators. 

Another interesting feature of the set up is that for small bare masses the pseudo-Goldstone 
modes satisfy a modified Gell-Mann-Oaks-Renner relation [12]. Indeed as it can be seen from 
the plot in figure 14 we have a linear relation (as opposed to square root one) between the 
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frequency uj and the bare mass m. It turns out that the slope of the observed hnear relation is 




Figure 14: Plots of the linear relation a) oc m and of the non-relativistic dispersion relation of the 
pseudo-Goldstone modes. 



given by: 



UJ = — Com , 

711] 



where cq represents the condensate at zero bare mass. In fact one can prove a more general 
result: 

LO = •yk -\ Cqui , (89) 

nr] 



where 



^ nri J ^(4f4 + 1)2 + i6r/2^4 





(90) 



and we have defined k = R'^k/b. Let us first verify the dispersion relation from equation (89) 



numerically. To do this end we consider a D5-brane embedding with a very small bare mass 
771 ~ 0.00062. Next we generate a plot of u versus k for fixed rj = 4.5. The resulting plot is 



presented in figure 14 The black fitting curve corresponds to the relation: 

^2 



UJ 



0.0004326 + 0.1010fc 



(91) 



On the other side the expression for 7 form equation (90) is 7 

0.00062 we obtain Ado/n 



0.1011 and by calculating 
0.0004321. One can see that the 



numerically Cq and for m 
relative error is 0.06% and 0.05% correspondingly. The analytic prove of the dispersion relation 



in equation (89) is a generalization of the zero temperature case considered in ref. [31]. Let us 
briefly provide the proof. 
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4.2.1 Low energy dispersion relation 



Let us consider the limit of small ui in equation ( 85 ) thus leaving only the linear potential term 
in uj. 



1 +/'2 



(92) 



Note that we have written down equation ( 92 ) in the dimensionless variables defined previously 
and used the notation P = P/i?^. It is convenient to define the following variables: 



1 + ['2 



(93) 



Then equation (92) can be written as: 



^ 0^ r^^+(4f4 + l)2 + l6r/2f4M 62-"- 



(94) 



Where the overdots represent derivatives with respect to p. Now if we take the limit m — ^ 
and A; ^ we have that u and obtain that: 



w=0 



e 



is a solution to equation (94). Our next step is to consider small m and expand: 

^ = 6 + 5^; = 6 + 59, 



(95) 



(96) 



where the variations 6C, and 6Q are vanishing in the m limit. Then, to leading order in m 
(keeping in mind that ~ m andlc^ ~ m) we obtain: 



^ e I ' 4f4 + 12 + i6r^2^4 I 



(97) 



Now we multiply equation (97) by 9 and integrate along p. The result is: 



(05^ - 650 - (0^0 - 0^0) -^V dpP{p) - -^-fk = 





(98) 



Using the definitions of 0, P{p) and ^ and requiring regularity at infinity for 77+, one can show 



that the first term in equation (98) vanishes and that: 



(656 - 656) ^ = g5m ; j dpP{p) = -n/A 





(99) 
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and hence using that 5fh = m we obtain equation (89) which we duphcate below: 

f 4 . . 
a; = 7A; H cqitl , 



(100) 



Further analysis of the pseudo goldstone spectrum would involve derivation of the effective 
chiral action along the lines of ref. [31]. We leave such studies for future investigations. 



5 Conclusion 



In this paper we studied flavoured holographic 1+2 dimensional Yang-Mills theory at finite 
temperature in an external magnetic field. One of the main results of our studies was the 
construction of the phase diagram of the theory in Section 3. It seems that the observed 
structure of the phase diagram can be understood based on rather general grounds. At a weak 
magnetic field the observed quadratic behaviour in figure |4] can be understood as representing 
the fact that the theory is describable by a Born-Infield like action and hence the free energy is 
an even function of the magnetic field. The square root behaviour at large magnetic fields on 
the other side represents competition of energy scales and thus is consequence of the freezing 
effect that the magnetic field has on the phase transition. Perhaps this could explain why the 
phase diagram in the coordinates used in figure |3] has the same shape as the phase diagram of 
the 1+3 dimensional case studied in refs. [251 [2l]. The similarity with the phase diagram of the 
Gross- Neveu model in 1+2 dimensions studied in refs. [13,117] could perhaps also be interpreted 
along these lines. It would be interesting to use alternative non-perturbative techniques to study 
the phase diagram of the defect field theory holographically dual to our set up and compare 
with the result obtained via the AdS/CFT correspondence. 

Interesting feature of the theory is the observation that the magnetic field and the temperature 
have the same effect on the theory for large bare masses. This is to be contrasted to the results 
for the D3/D7 set up. More precisely this suggests that one can tune the sign of the condensate 
at large bare masses. Furthermore, one has analytic control on that regime of the theory. It is 
worth looking for possible applications of this property of the system. 

Another interesting feature of the theory is the enhanced relative jump of the entropy at zero 
bare mass. It would be interesting to employ alternative methods and explore the jump of 
entropy. For example by counting how the number of microstates changes. It is also intriguing 
that the ratio of the entropy before and after the phase transition seems to be numerically very 
close to 27r suggesting that perhaps this ratio can be obtained in a closed form. 
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A possible direction for future studies is the derivation of the temperature dependence of the 
condensate of the theory at zero bare mass, using alternative field theory approach. In general 
such studies could verify the validity of the holographic set up. 

An interesting property of the D3/D5 set up is that the holographic renormalization of the 
D5-brane action does not require additional counter terms due to the external electric or mag- 
netic field. This enables one to define unambiguously the magnetization of the theory. This 
observation suggests that D5-brane probing of more complicated asymptotically AdSs geome- 
tries would still have the same property. It would also be interesting to further analyze the 
temperature dependence of the diamagnetic response in the deconfined phase and attempt to 
model that dependence based on the control that we have on the temperature dependence of 
other properties of the theory, such as the conductivity. 

The meson spectrum of the theory could also be analyzed in more details. For example a 
more complete study of the spectrum of quasi-normal modes including non-zero bare masses. 
Finally, studies of the effective chiral action along the lines of the studies performed in [3l] are 
of potential phenomenological interest. 
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